Preheating in New Inflation 

Mariel Desroche/ Gary N. Felder/ Jan M. Kratochvil,^ and Andrei Linde^ 

^Department of Physics, Clark Science Center, 
Smith College Northampton, MA 01063, USA 

^Department of Physics, Stanford University, Stanford, CA 94305-^060, USA 

\^ (Dated: February 1, 2008) 

o 

O ■ Abstract 
(N 

Jlj ■ During the last ten years a detailed investigation of preheating was performed for chaotic inflation 

. 

and for hybrid inflation. However, nonperturbative effects during reheating in the new inflation 
scenario remained practically unexplored. We do a full analysis of preheating in new inflation, using 
a combination of analytical and numerical methods. We find that the decay of the homogeneous 
component of the inflaton field and the resulting process of spontaneous symmetry breaking in the 
simplest models of new inflation usually occurs almost instantly: for the new inflation on the GUT 
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l/^ ' scale it takes only about 5 oscillations of the field distribution. The decay of the homogeneous 
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inflaton field is so efficient because of a combined effect of tachyonic preheating and parametric 
_ , resonance. At that stage, the homogeneous oscillating inflaton field decays into a collection of 

^ I waves of the inflaton field, with a typical wavelength of the order of the inverse inflaton mass. This 



X 



stage usually is followed by a long stage of decay of the inflaton field into other particles, which 
can be described by the perturbative approach to reheating after inflation. The resulting reheating 
temperature typically is rather low. 
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I. INTRODUCTION 



Inflationary cosmology solves a number of problems in the big bang model and is well 
supported by observational evidence. According to this theory the universe at the end of 
the inflationary period consisted almost entirely of the homogeneous inflaton field. After 
inflation this field decayed into inhomogeneous fluctuations and other forms of particles and 
fields, ultimately giving rise to the forms of matter that make up the universe today. 

An understanding of this decay period, known as reheating, can provide crucial links be- 
tween the inflationary epoch and the subsequent thermalized hot big bang era. The reheating 
period may have involved high energy phase transitions, symmetry breaking, baryogenesis, 
and other effects that could have observable signatures and that could give us insights into 
physics at energies beyond the reach of accelerator experiments. 

Early discussions of reheating were based on the assumption that the homogeneous infla- 
ton field decayed perturbatively as a collection of particles The perturbative mechanism 
typically requires thousands of oscillations of the inflaton field until it decays into usual ele- 
mentary particles. More recently, however, it was discovered that coherent field effects such 
as parametric resonance can lead to the decay of the homogeneous field much faster than 
would have been predicted by perturbative methods, within few dozen oscillations j2|. These 
coherent effects produce high energy, nonthermal fluctuations that could have significance 
for understanding developments in the early universe, such as baryogenesis. This early stage 
of rapid nonperturbative decay may be followed by a period of slower, perturbative effects, 
so the rapid early stage is called preheating. 

In it was found that another effect known as tachyonic preheating can lead to even 
faster decay than parametric resonance. This effect occurs whenever the homogeneous field 
rolls down a tachyonic {{(PV/ dcj)^) < 0) region of its potential. When that occurs a tachyonic, 
or spinodal instability leads to exponentially rapid growth of all long wavelength modes 
{k"^ < \{(PV/d(f)'^)\). This growth can often drain all of the energy from the homogeneous 
field within a single oscillation. 

We are now in a position to classify the dominant mechanisms by which the homogeneous 
inflaton field decays in different classes of inflationary models. The simplest of these models 
can be broken into three classes: small field, or new inflation models Q], large field, or 
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chaotic inflation models [5|, and multi- field, or hybrid models j6|. 

In simple chaotic inflation preheating is generally dominated by parametric resonance, 
although there are parameter ranges where this can not occur. ^ Tachyonic preheating does 
not occur in these models because the curvature of the potential is always positive. In 
it is shown that tachyonic preheating dominates the preheating phase in hybrid models of 
inflation. In this paper we explore preheating in new inflationary models. We find that for 
almost all realistic parameters tachyonic preheating works simultaneously with parametric 
resonance. The resulting effect is very strong, so that the homogeneous mode of the inflaton 
field typically decays within few oscillations. 

We should emphasize, however, that this stage of reheating is not the last one. The 
inhomogeneities of the scalar field can be very long-living in the new inflation scenario; their 
decay can be described by the perturbative methods of Ref. 

We do not pretend that the taxonomy given above exhausts the space of all possible 
inflationary theories. One can imagine arbitrarily complicated models that combine aspects 
of small and large field models or theories with non-canonical kinetic terms, gravitational 
effects, and so on. Nonetheless, we believe that by understanding preheating in these three 
classes of theories we will have a good understanding of all of the simplest possibilities, and 
in many cases the more complex models will simply involve combinations of the effects seen 
in these simpler cases. 

In Section O we give an analytic investigation of tachyonic preheating in the simplest 
model of new inflation. Section IIIll describes parametric resonance which occurs in this 
model. A full treatment, which would unify these two effects and take into account other 
important effects which may occur at the later stages of preheating, requires lattice cal- 



describe our lattice simulations and summarize their results, including the growth of fluctua- 
tions, the properties of the resulting spectra, and the formation of domain walls. Finally, we 
conclude with a summary of preheating in new inflation models and suggestions for further 
work. 

^ We refer to "simple" models as ones where the field simply oscillates around a minimum at = 0, such as 
V =^m?(p'^ . In non-oscillatory models such as exponential potentials the situation is more complicated. 
See 13 for a discussion of preheating in these models. 



culations (sf, which we perform using the program LATTICEEASY 
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II. TACHYONIC PREHEATING IN NEW INFLATION 



In this paper we are considering tachyonic preheating in models of new inflation, i.e. 
models where the inflaton field rolls from a potential maximum at = to a minimum at 
a symmetry breaking value (p = v. For definiteness we focus primarily on the original new 
inflation model 1^] based on the Coleman- Weinberg potential jl^ 



V{4>) = h4>' 



V 4 
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but the results we present should apply to any new inflation model (with one small exception 
noted below.) The potential energy takes its maximum value Vq = (l/16)Af^ at = and 
vanishes at the minima at = ±v. We are going to use units Mp = 1, even though sometimes 
we will write Mp explicitly. 

During inflation the field is near the maximum. Roughly speaking, inflation ends when 
^ Neglecting factors of order unity (including logarithmic factors) this occurs at a 
field value 

00 ~ — • (2) 

Mp ^ ' 

While this result is approximate, we have verified numerically that for a wide range of values 
of V this field value corresponds to a time when kinetic energy is still orders of magnitude 
smaller than potential energy, so we can safely start our lattice simulations at this time 
without fear of missing important effects. 

First we consider the homogeneous rolling field, neglecting backreaction. The energy lost 
to Hubble friction in one oscillation can be estimated as 

^ H^A(j). (3) 

Up to factors of order unity we can approximate H ^ y/Vo/Mp, (j) ^ y/Vo, and Acf) ^ v, so 

AV«V'„^ (4) 

Approximating the potential near the top as a negative quartic potential we can use this 
formula to find the field value 0i after the first oscillation 

-A,S + A,t^V„-^«A„'-^ (5) 
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Continued application of this argument leads to the conclusion that after n oscillations the 
field value will be roughly 

/ \ 1/4 

We wish now to estimate the growth of fluctuations due to tachyonic preheating. We 
consider preheating to be finished when 6(j) ~ v, at which point it is meaningless to talk 
about oscillations of the homogeneous field. To estimate this growth we consider the growth 
of a single mode (pk, which to first order obeys the equation of motion 

k + (p' + miff) (t>k = 0, (8) 

where To^ff — {(PV/dclP') is the effective mass of the field as a function of the value of the 
homogeneous field (0) and p is the physical (not comoving) momentum. Because miff < 
modes with < \mlff\ will experience exponential growth. 

To estimate the speed of this growth we can use the following trick. If we consider a 
homogeneous field we can define a variable D = (p and write an equation of motion for D 

Thus modes with <^ l"^e//l obey the same equation of motion as the derivative 0. We 
know the solution to the equation for (p from energy conservation, though. Treating the 

potential as an inverted quartic and neglecting Hubble friction, oc 0^. The growth of the 
fiuctuations 50 is determined by the modes with -C |^e//l' ^-Iso grows proportionally 
to 02. 

At the end of inflation 

s^^H^y^^^. (11) 

Thus after the held falls to the minimum at = v 
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Note that as the field rolls back up the potential the fluctuations will not lose this growth 
because the waves are decoherent. 

This result says that for v < \/A, which is typically of order 10^^, preheating will complete 
in a single oscillation. For larger values of v we must consider subsequent oscillations. In 
each one Scf) grows by an amount (v/0„)^, which is roughly (Mp/v)^/^, so after n oscillations 
the fluctuations will have grown to roughly 

50 ^ VXMp (^] . (13) 

The number of oscillations required for S(j) to reach v in this approximation is 

n^^^. (14) 

For A = 10^^^, which corresponds to COBE normalization of density perturbations, and 
V ~ 10~^Mp, which corresponds to the GUT scale, preheating will complete in roughly 5 
oscillations. Of course when 50 ~ v the perturbative approximation will break down, but 
preheating will be over. As we will see, the numerical investigation conflrms this simple 
estimate. 

Qualitatively these results depend only on the fact that the inflaton feels an inverted 
polynomial potential as it rolls from its maximum at = to the minimum at <!) — v. Thus 
the same type of behavior would result from an inverted quartic with or without logarithmic 
corrections, an inverted cubic potential, or nearly any other symmetry breaking potential.^ 

These estimates suggest that for any f <^ 1 tachyonic preheating will drain the energy 
from the homogeneous, oscillating, inflaton fleld within a few oscillations. We will see in 
the next section that there are important corrections to these rough estimates, and for v 
sufficiently large tachyonic preheating does not operate, however for v ^ 10"^ the behavior 
is quahtatively as described by these analytical estimates. 



^ These results do not hold for the case of an inverted quadratic potential. However, if you assume that the 
potential near the maximum is an inverted quadratic then is constant as the field rolls, so inflation 
cannot end in this regime. 
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III. PARAMETRIC RESONANCE 



In addition to the tachyonic preheating, there is also a parametric resonance in this 
model. The most efficient part of this process occurs far from the minimum of the effective 

H' ential, where it was expected to happen in the earlier works on reheating in new inflation 

Indeed, let us find the range of the values of the field 0, for which the adiabaticity 
condition is violated and a broad parametric resonance takes place. 

The adiabaticity condition is w < w"^. Here w'^ = k'^ + V" = k"^ + 3A0^ (In ^ + 1)- Instead 
of investigating the general case, let us look at the vicinity of the point 0*, where V" = 0, i.e. 
at the point where the low-momenta modes change their nature from normal to tachyonic. 
At that point w'^ = /c^, whereas w = and 0^ ~ O.OSAf ^. This gives a broad range of 
momenta k for which the adiabaticity condition is violated: 

k ^0.5 V\v = 0.5 m , (15) 

where m = \fXv is the mass of the inflaton field in the minimum of V{(j)) at = f . This 
means that the parametric resonance is powerful enough to excite all modes and create 
particles with momenta smaller than (the half of) the mass of the inflaton field at the 
minimum of its effective potential. These momenta, however, are much higher than the 
smallest momenta excited by tachyonic preheating. 

But this is also the range of momenta where the tachyonic amplification of fluctuations 
is operative. In fact, when the scalar field moves in the region with < 0*, its mass 
is tachyonic, so many of the modes produced due to strong nonadiabaticity near = 0* 
continue growing exponentially when the field moves in the region < 0*. 

As a result, parametric resonance and tachyonic growth amplify each other, which should 
lead to an even faster decay of the homogeneous scalar field, with the spectrum of produced 
fluctuations spanning a large range of momenta, from the inflationary Hubble constant 
H ~ VXvjf to the mass of the field in the minimum of its effective potential ~ \/~\v. The 
homogeneous scalar field should decay within just a few of its oscillations. In this sense, 
new inflation takes an intermediate position between hybrid inflation, where the decay is 
purely tachyonic and takes just a single oscillation, and the simplest versions of chaotic 
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inflation scenario, where the decay is due to parametric resonance and occurs within few 
dozen oscillations. In what follows we will describe the results of numerical simulations of 
reheating in new inflation. 

IV. NUMERICAL CALCULATIONS 

While the estimates in the previous section provide a useful general description of the 
process of tachyonic preheating in new inflation, a full analysis including backreaction and 
rescattering effects requires a lattice simulation. We performed a series of simulations of the 
Coleman- Weinberg model (P) using the LATTICEEASY program P|. We used A = 10~^^, 
which is fixed by the normalization of the cosmic microwave background, but we varied v as 
a free parameter. In all of the figures in this section time is measured in units of {VXvy\ 
wave number k is measured in units of VXv, and the inflaton field is measured in units of 

V. Note that ^/Xv is the mass of the inflaton field in the minimum of the effective potential; 
we use the gravitational Planck mass Mp = 1/ VG. 




The number of gridpoints required to meet both of these requirements was quite large, 
which restricted us to doing one and two dimensional simulations. We verified for each 
case that the essential features were unchanged in going from one to two dimensions, so we 
tentatively infer that the same basic results would hold in three dimensions as well. In the 
appendix we give all of the parameters and other details about the simulations. 

Figure m shows the evolution of the mean value of the inflaton field for v = O.lMp. The 
Hubble friction has a slightly greater effect on the field than our rough analytical estimates 
suggested. The result of this small correction is that after falling from the top of the potential 
the mean of the field never comes back up into the tachyonic region. Consequently tachyonic 
preheating and the parametric resonance described in the previous section does not occur for 
this high value of the symmetry breaking scale. Fluctuations of the field in this simulation 
never grew from their initial vacuum values. The decay of the homogeneous component 
of the inflaton field in this regime occurs very slowly. In the beginning, it may occur due 
io a very narrow and inefficient parametric resonance, as described in Section IV of Ref. 
Q. but then ve,v rap.d.v th,s completely disappears, and decay continues due to 

perturbative effects described in 
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FIG. 1: Mean value of the inflaton field ((/>) (zero mode) for v = 10~^Mp. A horizontal line 
indicates the field value below which the potential is tachyonic (negative mass squared). The 
amplitude of the field here and in other figures is given in units of v, whereas the time t is given in 
units m^^ = (^y/Xv^ , where m is the mass of the scalar field near the minimum of the effective 
potential. 

<0>^ , <(50^ > 




FIG. 2: Squared mean (c/))^ and variance (50^) of the inflaton field for v = 10~^Mp. A horizontal 
line indicates the field value below which the potential is tachyonic (negative mass squared). As 
one can clearly see form this figure, the oscillations of the homogeneous component of the scalar 
field (j) are completely damped out after the first 5 oscillations, whereas the variance remains much 
smaller than 1, in units of v. This means that the process of spontaneous symmetry breaking in 
this scenario completes within 5 oscillations. 

For smaller values of the parameter v the situation changes dramatically. Figure |21 shows 
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the mean and variance of the inflaton field for v = 10~^Mp. In this plot you can clearly 
see that the mean field oscillates through the tachyonic region repeatedly, leading to rapid 
exponential growth of fiuctuations. After about five or six such oscillations the zero mode 
is effectively destroyed. We did simulations of values of v ranging from lO^^Mp to lO^^Mp 
and found that the cutoff below which tachyonic preheating was efficient was somewhere 
in between 10~^Mp and 10~^Mp. Thus for nearly all realistic values of the parameters, 
tachyonic preheating may be expected to dominate the initial decay of the homogeneous 
inflaton fleld in models of the type we are considering. In the rest of this section we will 
focus on results for v = 10~^Mp as an illustrative example, but we found these results to be 
typical for parameters for which tachyonic preheating was possible. 




L k 
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FIG. 3: Spectra of occupation number as a function of momentum k at different times for 
V = IQ-^Mp. 

To explore the decay of the homogeneous inflaton in greater detail we can consider the 
occupation number Uk of different modes. Figure El shows the development of this spec- 
trum over time. At the flrst time shown the fleld has been slowly rolling near the top of 
the potential. Only very long wavelength modes have been excited by this point. In the 
next frame we see continued growth of long-wavelength modes due to tachyonic preheating. 
Tachyonic preheating can only amplify modes with k < |me//|. The effective mass of the 
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field varies from zero at tlie maximum to a minimum value of \meff\ ~ 0.5-\/Af, so modes 
with momenta above this cutoff were not excited. Even for modes below the cutoff, lower 
k modes were excited more strongly than higher k modes because the field spent more time 
in the regime where these low k modes were tachyonic. 

Following this initial stage we see the growth of higher k modes, leading to the devel- 
opment of a peak around k ^ \/Xv. The modes in this peak region are marginally able 
to be excited by tachyonic preheating. Moreover, tachyonic preheating would be expected 
to produce a monotonically decreasing spectrum like the one in the second frame of figure 
El rather than a peak such as we see in the third frame. This peak must therefore be the 
result of parametric resonance, as described in Section IIIII As a further check, we examined 
whether the formation of this peak was purely attributable to parametric resonance pro- 
duced by the oscillating zero mode, or whether it was also influenced by scattering from the 
already-produced long-wavelength modes. To this end we solved the coupled equations for 
the evolution of the zero mode (neglecting backreaction) and a single mode from the peak. 




FIG. 4: Growth of a fluctuation in the peak (k ~ 0.58). The lower curves show the evolution of the 
zero mode and the upper curves show the occupation number of this mode. The solid (black) 
curves show results from LATTICEEASY. The thin (red) curves show results from a Mathematica 
calculation in which the zero mode was evolved with no backreaction. These results confirm 
our expectations that the growth of the modes with k ~ 0.5 occur due to parametric resonance 
practically independently of the tachyon preheating. 

Figure |3 compares the growth of this mode in the lattice simulation to its growth when 
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coupled only to the zero mode. The two plots are nearly identical until the point when 
backreaction from the long- wavelength modes significantly affects the zero mode. Thus the 
formation of the peak in the third frame of figure El is purely a result of parametric resonance 
resulting from the oscillations of the zero mode. Note, however, that this plot also shows 
the need for lattice simulations; a linear analysis would not show the effects of backreaction 
and the point at which the decay of the homogeneous mode completes. 

Figure |31 shows a zoomed in view of the growth of the mode depicted in figure lU The 
growth of the occupation number occurs when the field is high on the potential hill and this 
growth alternates with periods of unchanging occupation number as the field moves through 
the minimum. This figure can be compared with the images showing broad parametric 
resonance in j^. In both cases the growth occurs when w^e// passes through zero, but in the 
case considered in that occurred around the potential minimum, whereas here it occurs 
on the side of the potential hill. 

1 



800 850 900^ 

FIG. 5: This figure is the same as figure ^ except showing only the lattice results and zoomed in 
to show features more clearly. As we see, the exponential growth of the occupation numbers occurs 
each time when the average value of the scalar field becomes much smaller than v. 

In summary, the fiuctuations of the field pass through several distinct stages. Initially 
modes with k ^ are excited by tachyonic preheating. Later modes with k ^ O.^VXv 

are excited by parametric resonance with the still oscillating zero mode, producing a peak 
in the spectrum. Later still this peak smooths out due to rescattering, the oscillations of 
the zero mode cease, the peak in the spectrum is smoothed out, and we end up with a fiat 
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spectrum in the IR followed by an exponential cutoff near k = v, as shown in the final frame 
of figure El This is similar to the spectrum that would be produced by parametric resonance 
in a large-field inflationary model, but the decay is significantly more rapid in this case. As 
a result, the total energy available to the high energy modes is very large. 

These processes are illustrated by a series of four panels in Fig. 6, which show a two- 
dimensional distribution of the inflaton scalar field at different stages of our calculations. 
The simulation was done for v = 10~^Mp on a 2048 x 2048 lattice, and the output was 
averaged over 8 lattice points for display in the figure. 

The first of the four panels in Fig. 6 shows the distribution of the field (f) at the time 
t = 740, at the very early stages of rolling of the field towards the minimum of its effective 
potential. At that time the average amplitude of the field was still very small (p ~ 0.00755. 
The interval of change shown in this panel is 0.0001. As we see, the amplitude of fluctu- 
ations of this field at that time is well below this level: The field looks almost completely 
homogeneous, apart for small quantum fluctuations which did not grow much at that stage. 

The second panel shows the distribution of the field (p at the time t = 830. The field still 
continues rolling down during the first oscillation, but the average amplitude of the field at 
that time has grown up to ~ 0.15395. The interval of change shown in this panel is also 
0.0001, as large as in the previous panel. We see that the long wavelength perturbations 
have already grown up significantly, and have magnitude 6(f) ~ 0.00005. 

The third panel shows the distribution of the field (p at the time t = 850. At that time the 
field already made one oscillation and started its way down during the second oscillation; 
compare with Fig. 2. Its average amplitude at that time is ~ 0.39785. The interval of 
change shown in this panel is 0.0002, i.e. two times greater than in the previous panels. As 
we can see, the amplitude of the long wavelength perturbations experienced an additional 
growth, but new perturbations with a much shorter wavelength have grown up on top of the 
long wavelength perturbations. 

Finally, the last panel in Fig. 6 shows the distribution of the field t = 1000. At that time 
the short wavelength fluctuations completely dominate the field distribution. The average 
value of the scalar field (0) no longer oscillates. Its value is close to 1, though it is slightly 
smaller than 1 because of the partial symmetry restoration due to the contribution of the 
fluctuations (0^) to the effective potential. The energy of the field is concentrated in the 
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waves of the field with momenta k ^1 produced by amphfication of quantum fluctuations 
during preheating. 









FIG. 6: Lattice images of a two-dimensional simulation with v = W~^Mp. In the beginning of the 
process, the scalar field (j) is practically homogeneous due to inflation, see the first panel. When it 
starts falling down, tachyonic instability generates long wavelength fluctuations of the scalar held, 
shown by the second panel. Then the short wavelength fluctuations are generated, and eventually 
the held distribution stabilizes near the minimum of the effective potential, which corresponds to 
~ 1 in this flgure. 
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If one would use perturbation theory to estimate the total time which it takes for the 
decay of the homogeneous component of the scalar field following one would find that 
this decay occurs only after 0(A^^) ~ 10^^ oscillations for A ~ 10^^^. It is quite amazing 
therefore that the nonperturbative effects lead to this decay within only 5 oscillations. The 
resulting field distribution after 5 oscillations, shown in the last panel of Fig. 6, is quite 
different from what one could expect on the basis of the perturbative approach to reheating 
or on the theory of parametric resonance in models of chaotic inflation. 

As we see from Fig. 3, the occupation numbers of ^-particles produced by the decaying 
infiaton field are exponentially large. This means that the initially homogeneous field 
decays into semi-classical waves of the field (j). Since the main contribution to the occupation 
numbers of produced particles is given by particles (waves) with very small momenta, one 
could be tempted to conclude that the main reason of the rapid decay of the homogeneous 
component of the scalar field is the tachyonic preheating. On the other hand, the phase 
volume of these modes is relatively small. We have verified that for v ~ 10^^ the main 
reason for the rapid decay of the homogeneous field in our computer simulations is the 
production of the particles with momenta k ~ 0.5 due to parametric resonance combined 
with the tachyonic effects, see Section IIIII We expect, however, that for v <^ 10~^, the 
main mechanism of the decay of the homogeneous mode will be related to the tachyonic 
preheating described in Section ITTl 



V. THE FINAL STAGE OF REHEATING: DECAY OF THE INFLATON TO 
OTHER FIELDS 

So far we have described analytically and numerically how tachyonic preheating and 
parametric resonance cause the homogeneous infiaton field to rapidly decay into its own 
fluctuations. This process leads to the decoherence of the infiaton field and to spontaneous 
symmetry breaking, i.e. to the disappearance of the amplitude of oscillations of its homo- 
geneous component. However, this is only the first part of reheating, which involves the 
decay of the infiaton into other degrees of freedom. In this section we discuss what happens 
when the infiaton is coupled to an additional field. We focus on the simplest case, a g'^cp'^x^ 
interaction with a second scalar field. Our basic conclusions, however, should be valid for 
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a range of possible models including ones where reheating produces vectors and spinors as 
well as scalars. 

Adding this field to the model gives us the potential 

1 



v{<p) = 



In^ 

V 4 



+ \9'<P\' + ^A.;^ (16) 



It is known that when fiuctuations of one field, here 0, are excited to exponentially large 
occupation numbers, that will lead to rapid production of particles of other fields to which 
it is coupled In this case, however, this is complicated by the fact that the symmetry 
breaking vev of the inflaton field also produces an effective mass for the additional field %, 
which tends to suppress its production. The mass of (f) in the minimum of its potential 
is given by = \f\v^ while the mass of x is given by = gcj) = gv. The parametric 
resonance near the minimum of V{(j)) is rather narrow; it occurs only when > 2m^. Thus 
the condition for x to be efficiently excited is g"^ <^ A. 

But this means that the decay to particles x occurs only if the decay rate is strongly 
suppressed by the small coupling constant g"^ ^ A ~ 10^^^. 

Figure [7| shows the number density of particles of the (p and x fields for this model. In 
both cases we take v = 10~^Mp, which we know produces efficient tachyonic preheating. In 
the case where g"^ = \ we see very little production of x particles. When g"^ = A/10, the 
X field is driven up to exponentially large occupation numbers some time after the field </> 
grows. We also tested g^ = lOOA and as expected found virtually no growth of the x field. 




FIG. 7: Growth of the occupation numbers of (j) and x iii the model (|16|) . In both plots the upper 
curve shows n<^ and the lower curve shows n^. 

The results of our calculations show that even in the cases when the nonperturbative 
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effects lead to the exponentially growing number of particles this number always remains 
exponentially smaller than the number of particles (j). This means that nonperturbative 
effects in this simple model lead to a rapid decay of the nearly homogeneous oscillating field 
into particles or waves of this field, but not to the final stage of reheating when all particles 
of the field decay to other particles. 

This suggests that the decay of the field eventually ends up by a perturbative stage 
when the field no longer oscillates coherently, and therefore one can use an approximation 
where each particle of the field decays independently, according to elementary theory of 
reheating developed in In this case the final temperature of reheating is given by 

Tj. ~ -\/r, where F = = |^ is the rate of decay — > XX i^i our model Ql; m = a/Aw . 

Simple estimates based on the decay rates for this model calculated in Q show that for 
= 0(A) the perturbative decay requires 0(A~^) oscillations, i.e. about 10^^ oscillations 

in the model of new inflation with A ~ 10^^^. The resulting temperature of reheating for 

V ~ 10-3, ^2 ^ A is ^ 10^ GeV. 

Whereas this scenario of the last stage of reheating seems rather general, there could 
be some important exceptions. For example, some fields interacting with the inflaton field 
may remain light even though their interaction with the inflaton field is very strong. This is 
possible if their mass is protected by some symmetry or if their mass is small in the minimum 
of the effective potential due to some cancellation mechanism. In this case the decay rate of 
the inflaton field into such fields can be quite large, the final stage of reheating may also be 
very rapid, and the resulting reheating temperature can be much higher. 



VI. CONCLUSIONS 

The first papers on reheating of the universe after inflation were devoted to investigation 
of this process in the new inflation scenario Q]. It was assumed that after inflation the 
nearly homogeneous scalar field oscillates for a very long time and eventually decays into 
particles of other fields in a process which can be described by a simple particle-by-particle 
decay of the scalar field 0. 

The subsequent investigation of reheating in chaotic inflation and in hybrid inflation 
have shown that reheating may occur much faster, due to nonperturbative effects such as 
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parametric resonance j2| and exponential growth of tachyonic modes j3|. In this paper we 
study nonperturbative effects during reheating in new inflation. 

We have found that for the simplest models of new inflation with a large amplitude of 
symmetry breaking, v > 10~^, the nonperturbative effects are relatively insignificant. After 
the first oscillation, the field (p never acquires the tachyonic mass. There can be a short 
stage of narrow parametric resonance, but it is inefficient, it shuts down very quickly, and 
the system enters the stage of perturbative particle-by-particle decay. A description of this 
regime can be found in Section IV of Ref. ^. 

On the other hand, in the versions of new inflation with v <^ 10~^, which includes the 
original new inflation model j^, the nonperturbative effects can be very powerful. As we 
have shown, a combination of tachyonic preheating and parametric resonance leads to a 
rapid dumping of energy of the oscillating inflaton field. For example, for v ~ 10^'^ the 
homogeneous oscillating inflaton field completely decays within 5 or 6 oscillations. For 
smaller values of v, the decay occurs even much faster. 

However, during this decay, the homogeneous oscillating inflaton field (j) decays to de- 
coherent waves of the field 0; its decay to the matter particles x af fhis stage typically is 
rather inefficient. Therefore the rapid nonperturbative stage of preheating is not the end of 
the story but only a prelude to a lengthy stage of perturbative decay, which eventually gives 
the reheating temperature that could be calculated by the methods developed more than 20 
years ago Q]. 

In this respect, the new inflationary scenario differs strongly from the simplest versions of 
chaotic inflation where the effect of parametric resonance and the subsequent violent stage 
of rescattering of produced particles makes thermalization much faster [12]. It differs even 
more from the hybrid inflation, where the process of decay of the zero mode of the inflaton 
takes just a single oscillation, and the subsequent stage of perturbative decay can occur 
relatively fast because the corresponding coupling constants there can be relatively large. 

The main reason of the slow decay of the inflaton field is related to the fact that the 
fields with which the inflaton field (p interacts strongly are getting heavy because of this 
interaction (this does not happen in the simplest versions of the chaotic inflation scenario). 
As a result, the inflaton field can only decay to the particles with which it almost does not 
interact This strongly suppresses the decay probability and the resulting reheating 
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temperature in the simplest models of new inflation. However, even in such models the 
reheating temperature can be sufficiently high for the subsequent stage of the low-scale 
baryogenesis. 



APPENDIX: PARAMETERS OF THE LATTICE CALCULATIONS 

The lattice simulation results shown in this paper were all calculated using LAT- 
TICEEASY a publicly available program for simulating interacting scalar fields in an 
expanding universe. The model (0) is encoded in the "colemanweinberg.h" model file on 
the LATTICEEASY website [14]. The website contains documentation on using LAT- 
TICEEASY and details on the algorithms employed. In this appendix we simply report 
the parameters used for these runs. These parameters, along with the Coleman- Weinberg 
model file, should be sufficient for anyone to reproduce all of the results reported here. 

The lattice spacing had to be chosen to be small enough to include modes with wave- 
lengths k > \fXv, the effective mass at the minimum of the potential, while the total 
size of the box had to be large enough to include significant numbers of modes with 



k < y A ln(t>/Mp)f^/Mp, the tachyonic mass at the end of infiation. This moment, 
= v'^/Mp, was taken as the start of the simulations. The initial velocity (p was set to 
zero. Note, however, that the program uses rescaled variables 

A _ „3/2^ (^ 

'rprogram y \ / 

tprogram \/~Xvt, (A. 2) 



SO that in program units the initial values of (p and (p are set to v and v'^ySn/S. (This latter 
result accounts for the derivative of the scale factor as well as (p.) 

Most of the results in Section HVl are for a run with v = 10^^ Mp. For this run we used a 
one-dimensional grid of 16, 384 gridpoints, a total box size L of 1500, and a time step of 0.01. 
We ran the simulation to a time t = 1000 but only began recording spectra at t = 800, the 
point when the field first started growing significantly. Various two-dimensional simulations 
were used to crosscheck many of the results in this section. The particular one shown in 
figure 6 was done on a 2048 x 2048 lattice using a time step of 0.05 (other parameters as just 
listed above). Note that all of these values are in program units. For example, L = 1500 in 
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program units corresponds to a physical box size of 1500/(v^i;) ^ 5 X W^M-\ 

The results in Section |3 used the same parameters as the v = 10~^Mp run described 
above, but with the addition of an extra field coupled to with coupling = Xoi = A/10. 

Figure n] shows results from a run with V = 10~^Mp. For this run we used 512 gridpoints, 
a total box size L = 40, a time step of 0.01, and a final time t = 100. 
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